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free approach is a factorization approach but does not require the coprime factorizations of the plant. We 
present two criteria for feedback stabilizability for MIMO systems in which transfer functions belong to the 
total rings of fractions of commutative rings. Both of them are generalizations of Sule's results in [SIAM 
J. Control Optim., 32-6, 1675-1695(1994)]. The first criterion is expressed in terms of modules generated 
from a causal plant and does not require the plant to be strictly causal. It shows that if the plant is stabilizable, 
the modules are projective. The other criterion is expressed in terms of ideals called generalized elementary 
factors. This gives the stabilizability of a causal plant in terms of the coprimeness of the generalized elementary 
factors. As an example, a discrete finite-time delay system is considered. 
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1. Introduction. In this paper we are concerned with the coordinate-free approach to 
control systems. This approach is a factorization approach but does not require the coprime 
factorizations of the plant. 

The factorization approach to control systems has the advantage that it embraces, within 
a single framework, numerous linear systems such as continuous -time as well as discrete- 
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time systems, lumped as well as distributed systems, 1-D as well as n-D systems, etc.[14]. 
This factorization approach was patterned after Desoer et al. [3] and Vidyasagar et al. [14]. In 
this approach, when problems such as feedback stabilization are studied, one can focus on 
the key aspects of the problem under study rather than be distracted by the special features 
of a particular class of linear systems. A transfer function of this approach is given as the 
ratio of two stable causal transfer functions and the set of stable causal transfer functions 
forms a commutative ring. For a long time, the theory of the factorization approach had been 
founded on the coprime factorizability of transfer matrices, which is satisfied in the case 
where the set of stable causal transfer functions is such a commutative ring as a Euclidean 
domain, a principal ideal, or a Bezout domain. However, Anantharam in [1] showed that 
there exist models in which some stabilizable plants do not have right-/left-coprime factor- 
izations. 

Recently, Shankar and Sule in [10] have presented a theory of feedback stabilization for 
single-input single-output (SISO) transfer functions having fractions over general integral 
domains. Moreover, Sule in [11, 12] has presented a theory of the feedback stabilization 
of strictly causal plants for multi-input multi-output (MIMO) transfer matrices, in which 
transfer functions belong to the total rings of fractions of commutative rings, with some re- 
strictions. Their approach to the control systems is called a "coordinate-free approach"([10, 
p. 15]) in the sense that they do not require the coprime factorizability of transfer matrices. 

The main contribution of this paper consists of providing two criteria for feedback stabi- 
lizability for MIMO systems in which transfer functions belong to the total rings of fractions 
of commutative rings: the first criterion is expressed in terms of modules ((ii) of Theo- 
rem 3.3) and the other in terms of ideals called generalized elementary factors ((iii) of Theo- 
rem 3.3). They are more general than Sule's results in the following sense: (i) our results do 
not require that plants are strictly causal; (ii) we do not employ the restriction of commutative 



rings. Further, we will not use the theory of algebraic geometry. 

The paper is organized as follows. In § 2, we give mathematical preliminaries, set up the 
feedback stabilization problem, present the previous results, and define the causality of the 
transfer functions. In § 3, we state our main results. As a preface to our main results, we 
also introduce there the notion of the generalized elementary factor of a plant. In § 4, we 
give intermediate results which we will utilize in the proof of the main theorem. In § 5, we 
prove our main theorem. In § 6, we discuss the causality of the stabilizing controllers. Also, 
in order to make the contents clear, we present examples concerning a discrete finite-time 
delay system in § 3, 4, 5 in series. 

2. Preliminaries. In the following we begin by introducing the notations of commuta- 
tive rings, matrices, and modules, commonly used in this paper. Then we give in order the 
formulation of the feedback stabilization problem, the previous results, and the causality of 
transfer functions. 

2.1. Notations. 

Commutative Rings. In this paper, we consider that any commutative ring has the iden- 
tity 1 different from zero. Let TZ denote a commutative ring. A zerodivisor in IZ is an 
element x for which there exists a nonzero y such that xy = 0. In particular, a zerodivisor x 
is said to be nilpotent if x" = for some positive integer n. The set of all nilpotent elements 
in TZ, which is an ideal, is called the nilradical of IZ. A nonzerodivisor in IZ is an element 
which is not a zerodivisor. The total ring of fractions of IZ is denoted by !F(1Z). 

The set of all prime ideals of IZ is called the prime spectrum of IZ and is denoted by 
Spec 7^. The prime spectrum of IZ is said to be irreducible as a topological space if every 
non-empty open set is dense in Spec IZ. 

We will consider that the set of stable causal transfer functions is a commutative ring, 
denoted by A. Further, we will use the following three kinds of rings of fractions: 



(i) The first one appears as the total ring of fractions of A, which is denoted by J-'{A) or 
simply by JF; that is, JF = {n/d \n,d E A, dis a nonzerodivisor}. This will be considered 
to be the set of all possible transfer functions. If the commutative ring A is an integral 
domain, JF becomes a field of fractions of A. However, if A is not an integral domain, 
then JF is not a field, because any zerodivisor of JF is not a unit. 

(ii) The second one is associated with the set of powers of a nonzero element of A. 
Suppose that / denotes a nonzero element of A. Given a set Sf = {1, /, /^, . . .}, which is 
a multiplicative subset of A, we denote by Af the ring of fractions of A with respect to the 
multiplicative subset Sj; that is, Aj = {n/d \n E A, d E Sj}. We point out two facts here: 
(a) In the case where / is nilpotent, Af becomes isomorphic to {0}. (b) In the case where / 
is a zerodivisor, even if the equality a/1 = b/l holds over Af with a,b E A, we cannot say 
in general that a = b over A; alternatively, a = b + z over A holds with some zerodivisor z 
of A such that zf^ = with a sufficiently large integer co. 

(iii) The last one is the total ring of fractions of Af, which is denoted by J^{Af); that 
is, J-'{Af) = {n/d \n,dE Af, dis a nonzerodivisor of Af}. If / is a nonzerodivisor of A, 
J^{Af) coincides with the total ring of fractions of A. Otherwise, they may not coincide. 

The reader is referred to Chapter 3 of [2] for the ring of fractions and to Chapter 1 of [2] 
for the prime spectrum. 

In the rest of the paper, we will use TZ as an unspecified commutative ring and mainly 
suppose that TZ denotes either A or Af. 

Matrices. Suppose that x and y denote sizes of matrices. 

The set of matrices over TZ of size x x y is denoted by TZ^^^. In particular, the set of 
square matrices over TZ of size x is denoted by {7Z)x. A square matrix is called singular 
over TZ if its determinant is a zerodivisor of TZ, and nonsingular otherwise. The identity 
and the zero matrices are denoted by E^ and O^xy, respectively, if the sizes are required. 



otherwise they are denoted simply by E and O. For a matrix A over TZ, the inverse matrix 
of A is denoted by A^^ provided that det{A) is a unit of !F{1Z). The ideal generated by IZ- 
linear combination of all minors of size m of a matrix A is denoted by Imiii^)- 

Matrices A and B over IZ are right-coprime over IZ if there exist matrices X and Y 
over 7^ such that XA + Fi? = E. Analogously, matrices A and B over 7^ are left-coprime 
over IZ if there exist matrices X and F over IZ such that AX + BY = E. Note that, in 
the sense of the above definition, even if two matrices have no common right- (^left-j divisors 
except invertible matrices, they may not be right- (^ left- jcoprime over 7Z. (For example, two 
matrices [zi ] and [Z2] of size 1x1 over the bivariate polynomial ring M.[zi, Z2] over the 
real field M are neither right- nor left-coprime over M.[zi, Z2] in our setting.) Further, a pair 
{N, D) of matrices N and D is said to be a right-coprime factorization of P over IZ if (i) 
the matrix D is nonsingular over IZ, (ii) P = ND^^ over J-'{TZ), and (iii) N and D are 
right-coprime over TZ. Also, a pair (N, D) of matrices N and D is said to be a left-coprime 
factorization of P over IZ if (i) D is nonsingular over IZ, (ii) P = D^^N over J-'(JZ), 
and (iii) iV and Z) are left-coprime over TZ. As we have seen, in the case where a matrix 
is potentially used to express left fractional form and/or left coprimeness, we usually attach 
a tilde '~' to a symbol; for example N, D for P = D-^N and Y, X for YN + XD = E. 



Modules. For a matrix A over TZ, we denote by Mr{A) [Mc{A)j the 7^-module gener- 
ated by rows (columns] of A. 

Suppose that A, B, A, B are matrices over TZ and X is a matrix over J^(TZ) such that 
X = AB-^ = B-^A with B and B being nonsingular. Then the 7^-module M^([ A* 5* ]*) 
(Mc([y4 -B])) is uniquely determined up to isomorphism with respect to any choice of 
fractions AB^^ yB^^Aj of X as shown in Lemma 2.1 below. Thus for a matrix X over JF(7^), 
we denote by Tx^n and 'Wx,ii the modules Mr{[A^ 5* ] ) and Mc{\A B]), respectively. 
If 7^ = ^, we write simply Tx and Wx for 7x ^ and Wx,a^ respectively. We will use, for 
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Fig. 2.1. Feedback system S. 



example, the notations Tp, Wp, Tc, and Wc for the matrices P and C over JF. 

An 7^-module M is called /ree if it has a basis, that is, a linearly independent system 
of generators. The rank of a free 7^-module M is equal to the cardinality of a basis of M, 
which is independent of the basis chosen. An 7^-module M is called ;?rq/ecrive if it is a direct 
summand of a free 7^-module, that is, there is a module A^ such that M © A^ is free. The 
reader is referred to Chapter 2 of [2] for the module theory. 

Lemma 2.1. Suppose that X is a matrix over J^(7V) and is expressed in the form of 
a fraction X = AB^^ = B^^A with some matrices A, B, A, B over TZ. Then the TZ-module 
Mr{[A^ B^] ) yMc{[A B])j is uniquely determined up to isomorphism with respect to 
any choice of fractions AB^^ yB^^Aj of X. 

Proof. Without loss of generality, it is sufficient to show that Mr{\A\ hiE] ) ~ 
Mr{\A\ B2] ), where bi e 71 and Ai{biE)^^ = ^2-82^^ Since hi is a nonzerodivisor 
andE2isnonsingular,wehaveM^([A* 61^]*) ^ M^([A* biEfB2)^Mr{[Al B\]^hiE) 
Mr {[A2 -82])- The other isomorphism can be proved analogously. D 

2.2. Feedback Stabilization Problem. The stabilization problem considered in this pa- 
per follows that of Sule in [11] who considers the feedback system S [13, Ch.5, Figure 5.1] 
as in Figure 2.1. For further details the reader is referred to [13]. Let a commutative ring A 
represent the set of stable causal transfer functions. The total ring of fractions of A, denoted 
by JF, consists of all possible transfer functions. The set of matrices of size x x y over A, 



denoted by A^^'^, coincides with the set of stable causal transfer matrices of size x x y. Also 
the set of matrices of size x x y over JF, denoted by J^^^v ^ coincides with all possible trans- 
fer matrices of size x x y. Throughout the paper, the plant we consider has m inputs and n 
outputs, and its transfer matrix, which itself is also called simply a plant, is denoted by P 
and belongs to jF"^™. We will occasionally consider matrices over A (jFJ as ones over Af 
or JF yJ^[Af)) by natural mapping. 
Definition 2.2. Define Fad by 

Fad = {{X, Y) e J^^"^^ X J^^^^ I det(F^ + XY) is a unit of T , 

X and y are positive integers). 

For P e J^"^™ and C G J^™^", the matrix H{P, C) E (J^)„+„ is defined by 
(2.1) H{P,C): 



{En + PC)-^ -P{Ern + CP)-' 

C{E^ + PC)-' {E^ + CP)-' 



provided {P, C) E Fs^d- This H{P,C) is the transfer matrix from [u\ u\] to [e\ e\] of 
the feedback system S. If(i) (P, C) G Fad <^nd (ii) H{P, C) G {A)m+n, then we say that the 
plant P is stabilizable, P is stabilized by C, and C is a stabilizing controller of P. 

2.3. Previous Results. In [11] Sule gave the results of the feedback stabilizability. We 
show them after introducing the notion of the elementary factor which is used to state his 
results. 

Definition 2.3. (Elementary Factors [11, p. 1689]) Assume that A is a unique fac- 
torization domain. Denote by T the matrix [ N* dE^ ] and by W the matrix [ N dEn ] 
over A, where P = Nd-^ with N G ^"^™, d e A Let {Ti, T2, . . . , TJ be the family of all 
nonsingular m x m submatrices of the matrix T, and for each index i, let fi be the radical of 
the least common multiple of all the denominators ofTT^^. The family F = {/i, /2, . . . , fr} 
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is called the family of elementary factors of the matrix T. Analogously let {Wi,W2, ■ ■ ■ ,Wr} 
be the family of all nonsingular n x n submatrices of the matrix W, and for each index j, 
let Qj be the radical of the least common multiple of all the denominators ofW^^W. Let 
G = {gi,g2, ■ ■ ■ ,gi} denote the family of elementary factors of the transposed matrix W*. 
Now let H = {fiQj \i = 1, . . . , r, j = 1, . . . , /}. This family H is called the elementary 
factor of the transfer matrix P. 

Then, Sule's two elegant results can be rewritten as follows. The first result assumes 
that the prime spectrum of A is irreducible. The second one assumes that ^ is a unique 
factorization domain. 

Theorem 2.4. (Theorem 1 of [11]) Suppose that the prime spectrum of A is irreducible. 
Further suppose that a plant P of size nxmis strictly causal, where the notion of the strictly 
causal is defined as in [12] (rather than [11]). Then the plant P is stabilizable if and only if 
the following conditions are satisfied: 

(i) The module Tp is projective of rank m. 

(ii) The module Wp is projective of rank n. D 

Recall that for a matrix X over JF, we use the notations Tx and Wx to denote ^-modules 
generated by using the matrix X. Further it should be noted that the definitions of Tp and 
yVp in [1 1] are slightly different from those of this paper. Nevertheless this is not a problem 
by virtue of Lemma 2.1. 

Theorem 2.5. (Theorem 4 of [11]) Suppose that A is a unique factorization domain. 
The plant P is stabilizable if and only if the elementary factors of P are coprime, that is, 
EneHih) =A. U 

2.4. Causality of Given Plants. Here we define the causality of transfer functions, 
which is an important physical constraint, used in this paper. We employ the definition 



of causality from Vidyasagar etal.\\A, Definition 3.1] and introduce two terminologies used 
later frequently. 

Definition 2.6. Let Z be a prime ideal of A, with Z ^ A, including all zerodivisors. 
Define the subsets V and Vs of T as follows: 

V = {n/d eJ^\neA, de A\Z], Vs = {n/d eJ^\neZ, de A\Z]. 

A transfer function in V \Psj is called causal (^strictly causalj. Similarly, if every entry of 
a transfer matrix over T is in V \Ps)y the transfer matrix is called causal (^strictly causalj. 
A transfer matrix is said to be 2-nonsingular if the determinant is in A\Z, and ^-singular 
otherwise. 

In [14], the ideal Z is not restricted to a prime ideal in general. On the other hand, 
in [11], the set of the denominators of causal transfer functions is a multiplicatively closed 
subset of A. This property is natural since the multiplication of two causal transfer functions 
should be considered as causal one. Note that this multiplicativity is equivalent to Z being 
prime provided that Z is an ideal. By following the multiplicativity, we consider in this paper 
that Z is prime. 

In this paper, we do not assume that the prime spectrum of A is irreducible and the 
plant P is strictly causal as in [11]. Alternatively, in the rest of the paper we assume only the 
following: 

Assumption 2.7. The given plant is causal in the sense of Definition 2.6. D 

One can represent a causal plant P in the form of fractions P = ND^^ = D~^N, where 
the matrices N, D, N, D are over A, and the matrices D, D are Z-nonsingular. 

3. Main Results. To state our results precisely we define the notion of generalized ele- 
mentary factors, which is a generalization of the elementary factors in Definition 2.3. Then 
the main theorem will be presented. 
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Generalized Elementary Factors. Originally, the elementary factors have been defined 
over unique factorization domains as in Definition 2.3. The authors have enlarged this con- 
cept for integral domains[8] and have presented a criterion for feedback stabilizability over 
integral domains. We enlarge this concept again in the case of commutative rings. 

Before stating the definition, we introduce several symbols used in the definition and 
widely in the rest of this paper. The symbol X denotes the set of all sets of m distinct integers 
between 1 and m + n (recall that m and n are the numbers of the inputs and the outputs, 
respectively). Normally, an element of X will be denoted by /, possibly with suffixes such 
as integers. We will use an element of X as a suffix as well as a set. For / G X, if ii, . . . , i^ 
are elements of / in ascending order, that is, «„ < 4 if a < h, then the symbol A/ denotes 
the matrix whose (A;, ?fc)-entry is 1 for i^ E I and zero otherwise. 

Definition 3.1 . (Generalized Elementary Factors) Let P G J^"-^"^, and N and D are 
matrices over A with P = ND~^. Denote by T the matrix [N^ D^] . For each / G X 
define the ideal Kpi of A by 

Api = {\eA\^Ke ^(™+")x'" AT = K/\iT}. 

We call the ideal Apj the generalized elementary factor of the plant P with respect to I. 

Whenever we use the symbol A with some suffix, it will denote a generalized elementary 
factor. We will also frequently use the symbols A and A/ with / G X as particular elements 
of Apj. Note that in Definitions 2.3 and 3.1, the matrices represented by T are different in 
general. However this difference is not a problem since the generalized elementary factors 
are independent of the choice of the fractions ND^^ as shown below. 

Lemma 3.2. For any I E T, the generalized elementary factor of the plant P with 
respect to I is independent of the choice of matrices N and D over A satisfying P = ND^^. 

Proof. Let A^, A^', D be matrices over A and d! be a scalar of A such that P = ND~^ = 
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N'd'~^ hold. Further, let 

In order to prove this lemma it is sufficient to show that the ideals Ap/i and Ap/2 are 
equal. Suppose that A is an element of Ap/i. Then there exists a matrix K such that 
A [A^* D*]* = K^i\N^ D^]\ Multiplying d'Em on the right of both sides, we have 
A [ A^'* d'Em ]^ D = KAj [ N'^ d'E„, ]* D. Since the matrix D is nonsingular, we have 
A [A^'* d'E^]^ = KAi[N^ d'Emf, so that A G Ap/2, which means that Ap/i C Ap/2. 
The opposite inclusion relation Ap/i D Ap/2 can be proved analogously. D 

Note also that for every / in X, the generalized elementary factor of the plant with respect 
to / is not empty since it contains at least zero. In the case where the set of stable causal 
transfer functions is a unique factorization domain, the generalized elementary factor of the 
plant with the matrix A/T being nonsingular becomes a principal ideal and the generator 
of its radical an elementary factor of the matrix T (in Definition 2.3) up to a unit multiple. 
Analogously, the elementary factor of the matrix W (in Definition 2.3) corresponds to the 
generalized elementary factor of the transposed plant P*. 

Main Results. We are now in position to state our main results. 

Theorem 3.3. Consider a causal plant P. Then the following statements are equiva- 
lent: 

(i) The plant P is stabilizable. 

(ii) The A-modules Tp and Wp are projective. 
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(iii) The set of all generalized elementary factors of P generates A; that is, 

(3.1) ^Ap, = A 

lei 

U 

In the theorem, (ii) and (iii) are criteria for feedback stabilizability. Comparing the the- 
orem above with Theorems 2.4 and 2.5, we observe the following: (ii) and (iii) can be con- 
sidered as generalizations of Theorems 2.4 and 2.5, respectively. For (ii), we do not assume 
as mentioned earlier that the prime spectrum of A is irreducible and the plant P is strictly 
causal. The rank conditions of Tp and Wp are deleted. For (iii), the commutative ring A is 
not restricted to a unique factorization domain. The elementary factors are replaced by the 
generalized elementary factors. Although two matrices T and W in Definition 2.3 are used 
to state Theorem 2.5, only one matrix T in Definition 3.1 is used in (iii). 

We will present the proof of the theorem in § 5. 

To make the notion of the generalized elementary factors familiar, we present here an 
example of the generalized elementary factors. 

Example 3.4. On some synchronous high-speed electronic circuits such as computer 
memory devices, they cannot often have nonzero small delays (for example [5]). We suppose 
here that the system cannot have the unit delay as a nonzero small delay. Further we suppose 
that the impulse response of a transfer function being stable is finitely terminated. Thus 
the set A becomes the set of polynomials generated by z^ and z^, that is, A = Mf^;^, z^], 
where z denotes the unit delay operator. Then A is not a unique factorization domain but 
a Noetherian domain. The total field JF of fractions of A is M^z"^ , z"^), which is equal to 'R(z) . 
The ideal Z used to define the causality is given as the set of polynomials in M[z^, z^] whose 
constant terms are zero; that is, Z = z'^A + z^A = {az"^ + bz^ \a,b E A}. Thus the set of 
causal transfer functions V is given as n/d, where n, d are in A and the constant term of d is 
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nonzero; that is, P = {n/iya + hz^ + cz'^) \n E A, a E R\{0}, b,cE A}. Further the set of 
strictly causal transfer functions Ps is given as Ps = {{aiz'^ + biZ^)/{a2 + b2z'^ + C2Z^) \a2 G 
M\{0}, ai,bi,b2,C2eA} 

Since some factorized polynomials are sometimes expressed more compactly and easier 
to understand than the expanded ones, we here introduce the following notation: a polyno- 
mial in M.[z] surrounded by "(" and ")" indicates that it is in A as well as in M.[z] even though 
some factors between "(" and ")" may not be in A. 

Let us consider the plant below: 



(3.2) 



P:-- 



(l-z3)/(l-z2) 
[l-8z^)/{l-Az^) 



eV 



2x1 



The representation of the plant is not unique. For example, the (l,l)-entry of the plant has 
an alternative form {I + z"^ + z'^)/{l + z^) different from the expression in (3.2). Consider 
parameterizing the representation of the plant. To do so we consider the plant P over M.(z) 
rather than over JF. Thus P can be expressed as 



(3.3) 



P 



{l + z + z^)/{l + z) 
'l + 2z + 4:z^)/{l + 2z) 



over KU 



However the coefficients of all numerators and denominators in (3.3) of z with degree 1 are 
not zero. To make them zero, we should multiply them by (ai(l — z) + biz'^ + ciz^) or 
(02(1 — 22;) + 62-2^ + C2Z^) with Oi, 61, Ci, 02, 62, C2 G ^ as follows 



(3.4) 



P 



{{l+z+z^){ai{l-z)+biz^+ciz'^)) 
((l+z)(ai(l-2)+fei22+ciz3)) 

{{l+2z+4z^){a2{l-2z)+b2z2+C2z'i)) 
{{l+2z){a2{l-2z)+b2z2+C2Z-i)) 



Every expression of the plant is given in the form of (3.4) with ai, 61, ci, 02, &2, C2 in A 
provided that the denominators are not zero. From this, we have two observations. One is 
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that the plant P does not have its right- and left-coprime factorizations over A (even so, it will 
be shown later that the plant is stabilizable). The other is that the elementary factor of this 
plant cannot be consistently defined over A. Thus we employ the notion of the generalized 
elementary factor. 



In the following, we calculate the generalized elementary factors of the plant. We choose 
the following matrices as N, D, and T used in Definition 3.1: 



111 



N :-- 



{1 - z'^){l - 4z^) 
'd]:=D:=[{l-z^){l-Az^)], T:=[N' D'f . 



n2 



Since m = 1 (the number of inputs) and n = 2 (the number of outputs), the set X is given as 
J = {{1}, {2}, {3}} and we let h = {1}, h = {2}, h = {3}. 

Let us calculate the generalized elementary factor Ap/^. Let ii = 1 so that Ji = {ii}. 
Then the (1, ii)-entry of the matrix A/^ is 1 and the other entries are zero. Thus we have 
A/^ = [ 1 ]. The generalized elementary factor Ap/^ is originally given as follows: 



(3.5) 



Apj^ = {\eA\3K e ^(™+")x™ AT = KAj^T} 

= {A G^|3A;i,A;2 G^ A[n2 df = ni[ki kif}. 



Consider (3.5) over R[z] instead of A. Then the matrix equation in the set of (3.5) can be 
expressed as 



A 



(3.6) 



(1 - z){l + z){l - 2z){l + 2z + 4^2)- 
{l-z){l + z){l-2z){l + 2z) 

{l-z){l-2z){l + 2z){l + z + z^) 



k2 
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The set of A's such that there exist /ci, /c2 E M.[z] satisfying (3.6) is given as {(1 + 2z){l + 
z + z'^)a I a e M.[z]}, denoted by Li. Then the intersection of Li and A is given as follows: 

(3.7) LinA={{il + 2z){l + z + z'^){a{l - 3z) + bz^ + cz^)) eA\a,b,ceA} 

This is equal to Ap/^ as shown below. First it is obvious that Li D A D Ap/^. For each 
(1 + 2z){l + z + z^)(a(l — 3^;) + bz"^ + cz^) with a,b,cE A, we have ki and k2 as follows 
from (3.6): 

ki = {l + z){l + 2z + 4z2)(a(l - 3z) + bz^ + cz^), 
k2 = {l + z){l + 2z){a{l - 3z) + bz^ + cz^). 

Both ki and k2 are in A. Hence LiCiAc Ap/^ and so Li fl ^ = Ap/^. By (3.7), we can also 
consider that Ap/^ is generated by ((1 + 2z){l + z + z^){l- 3z)), ((1 + 2z){l + z + z'^)z'^), 
dind{{l + 2z){l + z + z'^)z^). 



Analogously, we can calculate the generalized elementary factors Ap/j and Ap/3 of the 
plant with respect to I2 and Is as follows: 

Ap72 = {((1 + z){l + 2z + Az^){a{l - 3^) + bz^ + cz^)) eA\a,b,ceA}, 
Ap73 = {((1 + z){l + 2z){a{l - 3z) + bz^ + cz^)) eA\a,b,ceA}. 

Observe now that 

Ap,, 3 ((1 + 2z){l -3z){l + z + z^)) =: Ao/^, 
Kpi, 3 ((1 + z){l + 2z + 4z^){l -3z + z^)) =: Aq/^ 
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and further 



where 



C^h^Oh + Cih'^Oh — 1) 



4233-236462^ -39836z^ -201780z"^ - 113016z^ +75344z^ 



'^■fi " 5852 



G A, 



a/2 



10085+184182^+12114023+i3i852z"^+1130162S 
5852 



eA. 



Now let 
(3.8) 



A/i :— a/iAo/i G Ap/j, A/2 :— a 1^X012 G Ap/j. 



Thus Ap/^ + Ap/2 = A and A/^ + A/j = 1. Hence by Theorem 3.3, the plant P is stabilizable. 
D 

4. Intermediate Results. In this section we provide intermediate results which will be 
used in the proof of our main theorem stated in the preceding section. This section consists 
of three parts. We first show that a number of modules generated from plants, controllers, 
and feedback systems are isomorphic to one another. Next we develop the results which will 
help to show the existence of a well-defined stabilizing controller. We then give the coprime 
factorizability of the plant over Af, where / is an element of the generalized elementary 
factor of the plant. 

Relationship in terms of Modules between Transfer Matrices P, C, and H{P,C). 
The first intermediate result is the relations, expressed in terms of modules, among the ma- 
trices P, C, and H(P, C) as well as their transposed matrices. A number of modules are 
isomorphic to one another as follows. 

Proposition 4. 1 . Suppose that P and C are matrices over J^(JV). Suppose further that 
det{En + PC) is a unit of J^ilZ). 
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(i) The following IZ-modules are isomorphic to one another: 

'Tp,TZ © Ten, TH{P,C),n> T^H{P\C^^)\Tly y^H{P':,Ct),-R> '^H(C,P),Tl- 

(ii) The following TZ-modules are isomorphic to one another: 

Wp^Ti © Wen, Wh(p,c),'r., WH{pt,cty,'R., Th{p\c^)^'r., Wh(c,p),'r.- 
Further for a matrix X over J^{TZ), 

(iii) Tx^n ^ Wx\n and Wx,n ^ Txt^n- 

Note here that in the proposition above, the controller C need not be a stabilizing con- 
troller. For the case where C is a stabilizing controller, see Lemma 2 of [11]. 

We can consider that the proposition above, especially the relations Tp^-ji © Tc^n ~ 
'^H{p,c),n — '^H(c,p),n^ gives an interpretation of the structure of the feedback system in 
the sense that the module generated by the feedback system is given as the direct sum of the 
modules generated by the plant and the controller. In the proof ("(i)^(ii)") of Theorem 3.3, 
this proposition will play a key role. 

Proof. We first prove (iii). Let A and B be matrices over TZ with X = AB^^ . Then 
we have Tx,n - M^([A* fi*]*) ~ M^{\A* 5*]) ~ W[B-^yA\n - Wx\n- The other 
relation 'Wx,n — Txt,n can be proved in a similar way. 

Next we prove (i). Suppose that det(-E„ + PC) is a unit of J^(1Z). We prove the 
following relations in order: (a) Tp^n © ^c,7^ ^ Th(p,c),'R., (b) TH(P,c),n ^ ^H{p^c*)^7^, 
(c) Tuip^c^Y^n — yVH{p\c^),n^ (d) Th(p^c),'r. — '^h{c,p),'r.- 

(a) of (i). The proof of (a) follows mainly the proof of Lemma 2 in [11]. By virtue of 
Lemma 2.1, it is sufficient to show the relation Tp^n © ^,7^ — Mr{[Nj^ duEm+n] ) with 
Nh G {lZ)m+n, dff E TZ, where H(P, C) = Nffdjj^. Let A^, A^^ be matrices over TZ and d. 
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dc be in 7^ with P = Nd'^ and C = Ncd'^}. Further, let 



Q 



dcEn N 

-Nc dErr 



S 



dcEn O 
O dE„i 



From these we have Tpn^Tc^iz ~ Mr{[Q^ S^] ). In addition, since det{En+PC) is a unit 
ofJ^(7^), the matrix iV// is nonsingular, so that M^([Q* S^f) ~ M^([Q* S^f{det{NH)Em+n)) 
holds. A simple calculation shows that 



Q 

s 



idet{NH)E„ 



dnEm+n 
Nh 



adj(iV^)5. 



Because both matrices S and adj(A'^^) are nonsingular, we finally have that 



'^p,n®T'c,n 



Mr 



'Q' 


)^Mr{ 


'Q' 


.s _ 




.s _ 



{det{NH)Em+n)) ^ Mr 



dHErn+n 

Nh 



Tr 



H(p,c),n- 



(b) of (i). Observe that the following relation holds: 



(4.1) 



t r^t\t 



H{P\C 






Em 


HiP,C) 


' 


Er, 


Er. 







Em 






Let Nh and dn be a matrix over TZ and a scalar of TZ, respectively, with H(P, C) = N^dj/. 
Then (4.1) can be rewritten as follows: 



t r^t\t 






Em 


Nh{ 


' 


Em 


dn)-' 


En 


_ 




Er, 


_ 





H{P\C 



Hence, we have matrices A and B over IZ with H{P , C ) = AB such that 



A 






Em 


Nh, 


B = 


' 


E„ 


En 









En 






d 



H- 
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This gives the relation Th{p,c),ii - TH(p\c^Y,n- 

(c) of (i). This is directly obtained by applying (iii) to the matrix H{P*, C^y. 

(d) of (i). Between the matrices H{P, C) and H{C, P), the following relation holds: 



HiCP) 






— En 


H{P,C) 


Em 








O Em 

-En O 



Letting A''^ and du be a matrix over TZ and a scalar of TZ with H{P, C) = Nndj/ as in (b), 



we have matrices N'^ and D'^j over TZ such that 



K 



D' 



H 



O -En 



Em O 



o 



o 



O -En 



E„ 



O 



Nh 

duEm+n 



holds. Since H{C, P) = N'^D'^ ^ holds and the first matrix of the right-hand side of the 
equation above is invertible, the relation Th(p,c),'R — '^H{c,p),n holds. 
Finally, arguments similar to (i) prove (ii). D 

Before moving to the next, we prove an easy lemma useful to give results for the trans- 
posed plants. 

Lemma 4.2. A plant P is stabilizable if and only if its transposed plant P^ is. Moreover, 
in the case where the plant P is stabilizable, C is a stabilizing controller of P if and only 
ifC* is a stabilizing controller of the transposed plant P*. 

Proof. (Only If) Suppose that a plant P is stabilizable. Let C be a stabilizing controller 
of P. First, (P*, C*) is in P^d, since (P, C) G Pad and det(P„ + PC) = det(P™ + P*C*). 
From (4.1) in the proof of Proposition 4.1, if P(P,C) e (^)„+„, thenif(P*, C*) e {A)m+n- 

(If) Because (P*)* = P, the "If" part can be proved analogously. D 
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Z-nonsingularity of Transfer Matrices. In order to prove the stabilizability of the given 
causal plant, which will be necessary in the proof of the main theorem (Theorem 3.3), we 
should show the existence of the stabilizing controller. To do so, we will need to show that 
the denominator matrix of the stabilizing controller is nonsingular. The following result will 
help this matter. 

Lemma 4.3. Suppose that there exist matrices A, B, Ci, C2 over A such that the 
following square matrix is Z -nonsingular: 



(4.2) 



A Ci 
B C2 



where the matrix A is square and the matrices A and B have same number of columns. Then 
there exists a matrix R over A such that the matrix A + RB is Z -nonsingular. 

Before starting the proof, it is worth reviewing some easy facts about the prime ideal Z. 

Remark A. A. (i) If a is in A\Z and expressed as a = 6 + c with h,c E A, then at least 
one of h and c is in A\Z. (ii) If a is in A\Z and h is in Z, then the sum a + 6 is in A\Z. (iii) 
Every factor in A of an element of A\Z belongs to A\Z (that is, \ia,h E A and ah E A\Z, 
then a,b e A\Z). U 

They will be used in the proofs of Lemma 4.3 and Theorem 3.3. 

Proof of Lemma 4.3 . This proof mainly follows that of Lemma4.4.21 of [13]. 

If the matrix A itself is Z-nonsingular, then we can select the zero matrix as R. Hence 
we assume in the following that A is Z-singular. 

Since (4.2) is Z-nonsingular, there exists a full-size minor of [A* 5* ] in A\Z by 
Laplace's expansion of (4.2) and Remark 4. 4(i,iii). Let a be such a Z-nonsingular full-size 
minor of [A* B^] having as few rows from B as possible. 

We here construct a matrix R such that dei{A+RB) = ±a+z with z E Z. Since A is Z- 
singular, the full-size minor a must contain at least one row of B from the matrix [A* B*] . 
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Suppose that a is obtained by excluding the rows ii, . . . , i^ of A and including the rows 
ji, . . . , jk of B, where both of ii, . . . ,ik and ji, . . . , jk are in ascending order. Now define 
R = (^jj) by Ti-^j-^ = ■ ■ ■ = n^j^ = 1 and r^ = for all other i, j. Observe that det {A + RB) 
is expanded in terms of full-size minors of the matrices [E R] and [A* fi* ] from the 
factorization A + RB = [E R][A^ 5* ] by the Binet-Cauchy formula. Every minor of 
[E R] containing more than k columns of R is zero. By the method of choosing the rows 
from [A* B*] for the full-size minor a, every full-size minor of [A* 5* ] having less 
than k rows of B is in Z. There is only one nonzero minor of[E R] containing exactly k 
columns of R, which is obtained by excluding the columns ii, . . . ,ikOf the identity matrix E 
and including the columns ji , . . . , j^ of i?; it is equal to ± 1 . From the Binet-Cauchy formula, 
the corresponding minor of [A* B^] is a. As a result, det(A+RB) is given as a sum of ±a 
and elements in Z. By Remark 4.4(ii), the sum is in A\Z and so is det{A + RB) . The matrix 
A + RB is now Z-nonsingular. D 

Coprimeness and Generalized Elementary Factors. We present here that for each non- 
nilpotent element A of the generalized elementary factors, the plant has a right-coprime fac- 
torization over A\. This will be independent of the stabilizability of the plant. 

Lemma 4.5. (cf. Proposition 2.2 of [9]) Let Api be the generalized elementary factor 
of the plant P with respect to / G X and further let \/Api denote the radical of Kpi (as an 



ideal). Suppose that A is in y/Api but not nilpotent. Then, the Ax-module Tp^Ax is free of 
rank m. 

Proof. Let T, A^, D be the matrices over A as in Definition 3.1. Recall that Tp^Ax denotes 
the ^A-module generated by rows of the matrix T. By Definition 3.1, there exists a matrix K 
over A such that YT = KAjT holds for some positive integer r. Then we have a factoriza- 
tion of the matrix T over Ax as T = {\'~^K)(AjT), where all entries of the matrix X'^K 
are in Ax- In order to show that the ^A-module Tp^Ax i^ ^^^ of rank m, provided that A is not 
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nilpotent, it is sufficient to prove the following two facts: (i) The matrix A/T is nonsingular 
over Ax. (ii) There is a matrix X such that the matrix [ \~^K X ] is invertible over A\ and 
the matrix equation T = [ X~''K X ] [ (A/T)* O ]* holds. 

(i). Observe that the matrix D is nonsingular over Ax as well as over A. Since D = 
^{n+i,...,m+n}T = (X^"^ A{n+i,...,m+n}K)(AiT) holds (uotc that the suffix of the symbol A 
is an ordered set of m distinct integers between 1 and m + n as before Definition 3.1), the 
matrix A/T is also nonsingular over Ax provided that A is not nilpotent. 

(ii). Let ii , . . . , i„ be n distinct integers in ascending order between 1 and m+n excluding 
the integers in /. Then let X be the matrix whose (ik, A;) -entry is 1 for each i/. and zero 
otherwise. Then the determinant of [ \~^K X ] becomes ±1 since the matrix A~''A/i^ is 
the identity matrix of {Ax)m- ^ 

The lemma below will be used in the proof ("(ii)^(iii)") of the main theorem by letting 
IZ = Af, where / is an element of the generalized elementary factor of the plant but not 
nilpotent. 



Lemma 4.6. IfTZ-module Tp^i (VVp,7^j is free of rank m [n], there exist matrices A 
and B yA and Bj over IZ such that (A, B) is a right-coprime factorization \{A^ B) is a left- 
coprime factorization ) of the plant P (G jF(7^)"^'"j over TZ. 

Proof. This lemma is an analogy of the result given in the proof of Lemma 3 of [ 1 1 ] . See 
the proof of Lemma 3 of [ 1 1 ] . D 

Example A.l . (Continued) Here we continue Example 3.4. Let us follow Lemmas 4.5 
and 4.6 with the plant of (3.2). Let the notation be as in Example 3.4. 

First we proceed along the proof of Lemma 4.5. As an example, we pick Ji G X as / 
and A/j G Ap/^ as A in the proof of Lemma 4.5. Recall that for each A G Ap/, there exists 
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a matrix K such that AT = KAjT holds. In the case of A/^ G Ap/^, the matrix K is given as 



(4.3) 



K 



k3 



A 



/i 



ai,{{l + z){l-3z){l + 2z + 4z^)) 
ai,{{l + z){l + 2z){l-3z)) 



Thus we have the factorization T = (A-''ir)(A/,T): 



[l-z^){l-4z^ 
[I - 8z'^){l - z^ 
[1 - z^){l - Az^) ] 



\-^ai, ((1 + z){l- ?>z){l + 2z + 4^2)) 



[(l-r^)(l-4.2)], 



where r = 1 and A/^T = [ (1 — -2'^)(1 — 4^;^) ]. As shown in part (i) of the proof of 
Lemma4.5, A/^T = [ (1 — z^){l — Az"^) ] is nonsingular. 

The matrix X in part (ii) of the proof of Lemma4.5 is given as X = [[j ^ °]* by letting 
ii = 2 and ^2 = 3 according to Ji = {1}. We can see that the matrix 



(4.4) 



[\-^K X' 



0' 



\-'ai,{{l + z){l-?>z)(l + 2z + Az')) 1 



A7>/,((l + z)(l + 2z)(l-3z)) 



1 



is invertible. Therefore the A\j -module Tpj^^ is free and its rank is 1. (However we will 
see that the .A-module Tp is not free, see Example 5.2) 



From (4.4) and the matrix equation T = [ X'^'K X ] [ {Ai^Tf O ]\ we let 



(4.5) 



N, 



h 



Du 



A7/aj,((l + z){l - 32)(1 + 2z + Az^)) 
A7V((l + ^)(l + 2^)(l-3^)) 



A7.'i^, 
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(4.6) 



X X 



1 

-A7>/,((1 + 2)(1-32)(1 + 2^ + 4^2)^ 1 



-A7/«7,((l + ^)(l + 2^)(l-3z)) 



1 



where Nj^ E A^^ , Yj^ e A^^ , Di^,Xi^ e {Axj )i, and x denotes some matrix. Then 
{Nj-^ , Di-^ ) is a right-coprime factorization of the plant over Axj^ with Y/^ Nj^ +Xi-^ Di^ = Ei . 
which is consistent with Lemma 4.6. D 

5. Proof of Main Results. Now we give the proof of the main theorem. 
Proof of Theorem 3.3. We prove the following relations in order: (a) "(i)^(ii)," (b) 
"(ii)^(iii)," and (c) "(iii)^(i)." 

(a) "(i)^(ii)": Suppose that C is a stabilizing controller of the plant P. Then, the ^-module 
Th(p^c) is obviously free. By the relation Tp^ji © Tc^-r ~ Th(p^c),ti in Proposition 4. l(i), we 
have that the ^-module Tp is projective. By using Proposition 4. l(iii) and Lemma 4.2, the 
projectivity of the ^-module Wp can be proved analogously. 

(b) "(ii)^(iii)": Suppose that (ii) holds, that is, the modules Tp and Wp are projective. We 
let T, A^, D be the matrices over A as in Definition 3.1. According to Theorem IV.32 of [6, 
p.295], there exist finite sets Fi and F2 such that (1) each of them generates A and (2) for 
any f E Fi [f e F2), the ^/-module 7p,^^ {y^p,Af) is free- Let F be the set of all /1/2's 
with /i G Fi and /2 G F2. Then F generates A again, and the ^j-modules Tpj^^. and 
'y^p,A[ are free for every f E F. We suppose without loss of generality that the sets Fi, F2, 
and F do not contain any nilpotent element because 1 + x is a unit of A for any nilpotent x 
(cf. [2, p. 10]). (However, we note that other zerodivisors cannot be excluded from the set F.) 
The rank of the free ^j-module Tp^Af is 'm, since m rows of the denominator matrix D are 
independent over Af as well as over A. Analogously the rank of Wp^Af is n. 

In order to show that (iii) holds, it suffices to show that the relation X]feF(/^) *^ 
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Yliiei^Pi holds for a sufficiently large integer ^. Once this relation is obtained, since 
Y^fepU^) = vA holds, we have Y.i&x^pi = ^■ 

Let / be an arbitrary but fixed element of F. Let Vf be a square matrix of size m whose 
rows are m distinct generators of the ^/-module Mr{[N^ /}*] ) (~ Tp^Aj)- We assume 
without loss of generality that Vf is over A, that is, the denominators of all entries of Vf 
are 1. Otherwise if Vf is over Af but not over A, Vf multiplied by /^, with a sufficiently 
large integer x, will be over A, so that we can consider such Vff^ as "V/." Thus the following 
matrix equation holds over A: 

(5.1) rT = KfVf 

with a nonnegative integer v and a matrix Kf E j[(m+n)xm_ 

In order to prove the relation ^ f^pif^) '^ Tliei ^^^' ^^ ^^^^ ^^^^ show the relation 

(5.2) /™^(r^/)c J]Ap, 

/ex 

and then 

(5.3) (/«) C IrnAirKf). 

Observe first that det{f''AjKf) E Apj because 

det{rAjKf)T = r''Kf8.dj{AjKf)AjT. 

Since every element of ImAif^f) i^ ^^ ^-linear combination of det{f''AjKfys for all 
/ G X, we have (5.2). 

We next present (5.3). Let A^o and Dq be matrices with Kf = [Nq Dq] ■ Since each row 
of Vf is generated by rows of [N^ D^] over Af, there exist matrices Yq and Xq over Af 
such that Vf = [ Y^f X^f ][N^ -D* ]*. Thus, since Vf is nonsingular over Af, we have 
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[Yq Xq][Nq D^] = Em, which implies that {Nq,Do) is a right-coprime factorization 
of the plant P over Af. Recall here that Wp^Af is free of rank n. Thus by Lemma 4.6 there 
exist matrices Nq and Do such that (Nq, Dq) is a left-coprime factorization of the plant P 
over Af. Let Yq and Xq be matrices over Aj such that NqYq + DqXq = En holds. Then we 
have the following matrix equation: 



(5.4) 



^0 -^0 

-Do No 



No -Xo 
Do Yo 



Em —^0-^0 + -^0^0 



O 



En 



Denote by R the matrix [ — Xq Yq] . Then the matrix [Kj R]is invertible over Af since 
the right-hand side of (5.4) is invertible. For each I E I, let / be the ordered set of n distinct 
integers between 1 and m + n excluding m integers in / and let ii, . . . , i„ be elements of / in 
ascending order. Let Aj G ^'"x(™+") denote the matrix whose (A;, ifc)-entry is 1 for ik E I 
and zero otherwise. Then, by using Laplace's expansion, the following holds: 



det{[Kf R]) = ^(±det(A/fs:/)det(Aji?)), 



lei 



which is a unit of Af. From this and since the ideal ImAf (Kf) is generated by det(A/i^j)'s 
for all Iel,we have ImAfi^f) = -^f- This equality over Af gives (5.3) for a sufficiently 
large integer ^. 

From (5.2) and (5.3), the relation J^fepif^) "^ J2i£i^Pi holds. Therefore we conclude 
that the relation Xl/ex ^pi ~ -^ holds. 

(c) "(iii)^(i)": To prove the stabilizability, we will construct a stabilizing controller of the 
causal plant P from right-coprime factorizations over Af with some /'s in A. Let N and D 
be matrices over A such that P = ND~^ and D is Z-nonsingular. From (3.1), there ex- 
ist A/'s such that ^ A/ = 1, where A/ is an element of generalized elementary factor Ap/ of 
the plant P with respect to / in X; that is. A/ G Ap/. Now let these A/'s be fixed. Further, 
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let X" be the set of J's of these nonzero A/'s; that is, Xl/extt ^^ — 1- ^^ ^^^ (t")' ^^ '-^'^ 
consider without loss of generality that for every / G X**, A/ is not a nilpotent element of A. 
For each / G X**, the ^A/-module Tp^Ax is free of rank m by Lemma4.5. As in (b) again, 
let Vxj be a square matrix of size m whose rows are m distinct generators of the ^Aj -module 
Mr{[N^ -D* ] ) (~ '^p,Ax ) ^^^ assume without loss of generality that Vxj is over A. Then 
there exist matrices X/, Y}, Nj, Dj over A\; such that 



YjNj + XjDj = E^ 



with P = NiDj^ over J^{Ax,). 

We here present a formula of a stabilizing controller which is constructed from the matri- 
ces Xj and Yj. For any positive integer cu, there are coefficients a/'s in A with Xl/extt ^t/^/ = 
1. Let cj be a sufficiently large integer. Thus the matrices X^DiXj and XfDiYi are over A 
for every I eXK The stabilizing controller we will construct has the form 



(5.6) C={J2 aiX'iDiXi)-\Y, ajX'^DjYj). 

/GJtt /GJtt 



In the following we first consider that the matrix (Xl/extt o/Aj-D/X/) is Z-nonsingular 
and show that the plant is stabilized by the matrix C of (5.6). After showing it, we will be 
concerned with the case where the matrix (X]/ejtt cljX'^DiXj) is Z-singular. 

Suppose that the matrix (X]/ejtt aiX^DjXi) is Z-nonsingular. To prove that C is a stabi- 
lizing controller of P, it is sufficient to show that (P, C) G Pad and that four blocks of (2.1) 
are over A. 
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We first show that (P, C) E Fad- The following matrix equation holds: 

E^ + CP = E^ + {J2 ajX^DjXj)-\J2 aiX':/DjYj)ND-' 
leii lexi 

= {J2 diXIDiXi)-^ 

7gj8 /extt 

By the (l,l)-block of (5.8), we have 

(5.7) E^ + CP= (J] ai\^,DiXi)-\ 

/gjti 

This shows that det(-E'm + CP) is a unit of T so that (P, C) G Pad- 

Next we show that four blocks of (2.1) are over A. The (2, 2)-block is the inverse of 

(5.7): 

{Ern + CP)-^ = Y,cii>HDiXi. 

/GX« 

Similarly, simple calculations show that other blocks are also over A as follows: 
(2, l)-block: 

C{En + PC)-^ = Y,ai\^,DiYi, 

/extt 



;i, l)-block: 



(P„ + PC)-^ = Er,-J2 ai>^1NiYi, 



I&V 



;i,2)-block: 



-P(P„ + CP)-^ = -J2 aiX'iNiXi. 

lexi 



To finish the proof, we proceed to deal with the case where the matrix (X]/ei8 a/Aj D/X/) 
is Z-singular. To make the matrix Z-nonsingular, we reconstruct the matrices X/ and Yi with 
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an 



lelK 



Since the sum of a/Ay s for / G XMs equal to 1, by Remark 4.4(i, iii) there exists at least 
one summand, say a/„A/ with an Jq G X", such that both a/^ and A/^ belong to A\Z. Let i?/^ 
be a parameter matrix of A^^^\ Then the following matrix equation holds over A\j : 



where Dj^ = det(A/jj-D/J-E„ and Nj^^ = X'^gNjg adj(Ayj|L'/J. Since a; is a sufficiently large 
integer, the following matrix equation is over A: 

{xI{Yj, + RjM){xiNj,) 

+{xi{Xj, - Rj^,N,,)){xiDj,) = Xj-E^, 



where the matrices surrounded by "(^" and "j" in the left-hand side are over A. From the 
first matrix equation of (5.5), det(Z}) = det(Dj) det{Vxj) over Axj for every / G X". 
Thus by Remark 4. 4(iii) the matrix A/y-D/j, is Z-nonsingular and so is the matrix D/(, (= 
detiXlDj,,)E^). 

Consider now the following matrix equation over A: 



(5.8) 



Z]/ejtt aiX'^DiXi 



E,ex« aiX^DiYr 



-ai.Xl det(A^ D,JiV,„ a,„A^, dei{XlDi,)Di, 

D 



D O 

N E„, 

[O aj,Xldet{XlDj,)Dj^ 



The (1, l)-block of (5.8) can be understood in the following way. From the last matrix 
equation in (5.5) we have the following matrix equation over Ax/. 



DjYjN + DjXiD = D. 
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Considering the above equation multiplied by a/Aj over A, we have the following equation 
over A: 

(5.9) aiX^jDiYiN + aiX'^DiXjD = ajX'^D + ajX'^Z, 

where Z is a matrix over A such that X^Z is the zero matrix for some positive integer x. 
Since c<j is a large positive integer, we can consider that the matrix ajXfZ in (5.9) becomes 
the zero matrix. Therefore, the (1, l)-block of (5.8) holds. Then the matrix of the right-hand 
side of (5.8) is Z-nonsingular since both of the matrices D and a/„A'^j det(Aj|^-D/,))i^/o in 
the right-hand side of (5.8) are Z-nonsingular. Hence the first matrix of (5.8) is also Z- 
nonsingular by Remark 4.4(iii). By Lemma 4.3 and (5.8), there exists a matrix R'j^ of A"^^" 
such that the following matrix is Z-nonsingular: 

(5.10) Yl (^i^'/DjXi - ar,XJ';;Dj^ a.dj{XlDjM^Nr,. 

Let now i?,„ := A^^, adj(Ay^ D,Ji?;,,, Xj, := Xj, - Ri.Ni,, and Yi, := Yj, + Ri.Dj,. Then 
the matrix J2iei» o-i^fDiXi becomes equal to (5.10) and Z-nonsingular D 

Remark 5.1. From the proof above, if (i) we can check (3.1) and if (ii) we can construct 
the right-coprime factorizations of the given causal plant over Axj for every I E X'^, then we 
can construct stabilizing controllers of the plant, where A/ is an element of the generalized 
elementary factor of the plant. For (i), if we can compute, for example, the Grobner basis [4] 
over A and if the generalized elementary factors of the plant are finitely generated, (3.1) 
can be checked. For (ii), it is already known by Lemmas 4.6 and 4.5 that there exist the 
right-coprime factorizations of the plant over Axj . 

Let us give an example concerning the Grobner basis. Consider the generalized elemen- 
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tary factors of Example 3.4. They are expressed as 

Api, = (((1 + 2z){l + z + z^){l- 3z))) + (((1 + 2z){l + Z + z^)z')) 

+ (((l + 2^)(l + z + ^V)), 
AP7, = (((1 + z)(l + 2^ + 4z')(l - 3^))) + (((1 + 2)(1 + 2^ + 4^2)^2^^ 

+ (((l + 2)(l + 22 + 42=')z=')), 

Ap,3 = (((1 + 2;)(1 + 2z)(\ - 3z))) + (((1 + z){l + 2z)z')) 

+ {{{l + z){l + 2z)z')). 

Hence each of them has three generators and so is finitely generated. Suppose here that 
we can calculate the Grobner basis over A (of Example 3.4). Then as above the plant is 
stabilizable if and only if the Grobner basis of the nine generators contains 1. D 

In the following two examples we follow the proof of Theorem 3.3. In the first one, we 
construct a stabilizing controller with part (c). In the other example, we follow part (b). On 
the other hand we do not follow part (a) since it can be followed easily with part (a) of (i) in 
the proof of Proposition 4. 1 . 

Example 5.2. (Continued) We continue Example 3.4 (and 4.7) and construct a stabiliz- 
ing controller of the plant as in "(iii)^(i)" of the proof above. Let the notation be as in 
Examples 3.4 and 4.7. 

Since, in this example, Ap/^ + Ap/j = A holds, X^ = {/i, I2}. For Ji G X", the matri- 
ces Nj-^, Dj^, Xj-^ and F/j of (5.5) over Axj have been calculated as (4.5) and (4.6). For 
I2 G X", the matrices Nf^, Dj^, Xf^ and Yi^ of (5.5) over Axj can be calculated analogously 
as follows: 



iV, 



h 



A7>7,((1 + 2z){l -3z + z^){l + z + z^)) 
1 
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Dj, = [ A7>,,((1 + z){l + 2z){l -3z + z-")) ] , 

17, = [0 i],x,, = [0]. 

Then the following matrices are over A: 

\i,Dj,Xj, = [0], Xj,Dj,Yj, = [aj,{il + z)il + 2z){l-3z)) 0], 
Xi,Dj,Xj, = [0], Xj,Dj,Yj, = [0 ai,{{l + z){l + 2z){l-3z + z^))]. 

Hence in this example, we can let c<j = 1 as a sufficiently large integer and a/ = 1 for all 
/ G X» (since E/ei« A^ = 1). 

Note here that the matrix Xj^Di^Xi^ + Xi^Di^Xj^ is Z-singular. Hence we should re- 
construct the matrices Y}- and Xj. with i being either 1 or 2 as in the proof of Theorem 3.3. 
Since, in this example, both A/^ and A/, are nonzerodivisors, we can choose each of 1 and 2. 
This example proceeds by reconstructing the matrices F/^ and X/^ , which means that Ji is 
used as Iq in the proof of Theorem 3.3. The actual reconstruction is done by following the 
proof of Lemma 4.3. 

Consider the first matrix of (5.8). Recall that iV/^ = A^^iV/o adj(A';'gL'/o) and 5/„ = 
det(AjgD/„)-E„. In this example, they are given as 



Ni, = {N, 



h 



_ai,{{l + z){l-3z){l + 2z + 4:z'')) 
bi, = {bi, =)ai,{{l + z){l + 2z){l - 3z))E2. 



One can check that the first matrix of (5.8) is Z-nonsingular. Then we construct a matrix R'^^^ 
of A^^"^ such that (5.10) is Z-nonsingular. To do so, we follow temporarily the proof of the 
Lemma 4.3. 
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Consider the first matrix of (5.8) as the matrix of (4.2), that is, 

B = -ai,\ldet{\lD,,)Ni, 

\ K 

= -ai,\i,{{l + z){l + 2z){l-^z)) 

_ai,{{l + z){l-?,z){l + 2z + Az^)) 

Then we choose a full-size a minor of [A* B* ] having as few rows from B as possible. In 
this example, we can choose both entries in B. Here we choose the (1, l)-entry of B, so that 

(5.11) a = -ai,\\{{l + z){l + 2z){l - 3z)). 

Thus we have A; = 1, ii = 1 and ji = 1, where the notations k,ii, . . . ,ik, and j'l, . . • , jk are 
as in the proof of Lemma 4.3 . Hence R in the proof is given as i? = [ 1 ] . We can confirm 
that A + RB = [a] which is Z-nonsingular by observing that every factor of the right-hand 
side of (5.1 1) has a nonzero constant term. 

From here on we proceed with following again the proof of Theorem 3.3. The notation R 
used above corresponds to the notation R'^^^ in the proof of Theorem 3.3 (that is, R'j^ = 
[1 ]). The matrix Rj^ is given as follows: 

Ri, = {Rj, =)\^^ adj(Ar,/^/Ji?^, = A,, [ 1 0] . 

Then new X/^ and F/^ are given as follows: 

Xi^ := Xj^ - Ri^Ni^ = [ -Ajj ] , 
Yj, := Yj, + Rj.Dj, = [ 1 + a,, A,, ((1 + z){l + 2z){l - 3z)) 0] . 
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Therefore a stabilizing controller C of the form (5.6) is obtained as 



C = {Xi,Dj,Xj, + Xj,Dj,Xj,)-\Xj,Dj,Yj, + Xi,Dj,Yj,] 
-1 



ai,Xl{{l + z){l + 2z){l-3z)) 



«/,((! + z){l + 2z){l - 3z)){l + ai,Xi,{{l + z){l + 2z){l - 3z))) " * 
ai,{{l + z){l + 2z){l-3z + z^)) 



The matrix H{P, C) with the stabilizing controller C above over A is expressed as fol- 
lows: 



H{P,C) 



hii hi2 his 

/l21 /i22 ^23 
hsi hs2 ^33 . 



where 



hn = -ai.Xliil + z){l + 2z){l - 3z)) 

+«/,((! + z){l -3z + z^){l + 2z + Az^)), 
hu = -ai,{{l + 2z){l + Z + z^){l -3z + z^)), 
hi3 = Xj^, 
h2i = -a/i((l + z){l - 3z){l + 2z + 4z^)) 

{l + Xi,ai,{il + z)il + 2z)il-3z))), 
h22 = «/,(((! + 2z){l - 3z){l + Z + z2))(l + a/,A/,((l + z){l + 2z){l - 3z))) 

~Xl{{l + z){l + 2z){l-3z))), 
h23 = ai.Xliil + z){l- 3z){l + 2z + Az^)), 

hi = aj,{{l + z){l + 2z){l - 3z)){l + aj,Xi,{{l + z)(l + 2z)(l - 3z))), 
h32 = «/2((l + ^)(1 + 2^)(1 - 3^ + z^)), 



35 

hss = -ai.Xliil + z){l + 2z){l - 3z)). 

Before finishing this example, let us show that the w4-module Tp is not free. We show it 
by contradiction. Suppose that Tp is free. Then the ^-module Mr{T) is also free. Since the 
matrix D, a part of T, is nonsingular, the rank of Mr{T) is m. Let V be a matrix in {A)m 
whose rows are m distinct generators of Mj.{T). As in (5.5), we have matrices Y , X, N', D' 
over A such that 

[iV* £)*]* = [A^'* £)'*]* \/, [Y X][N^ £)*]* = y, Y' N' + X' D' = El. 

However the last matrix equation is inconsistent with the fact that the plant P does not have 
coprime factorization. Therefore Tp is not free. Nevertheless we note that Tp is projective 
by Theorem 3.3. D 

Example 5.3. (Continued) Let us follow part (b) in the proof of Theorem 3.3. Suppose 
that (i) of Theorem 3.3 holds, that is, the modules Tp and Wp are projective. 

Consider again the plant P of (3.2). Let Fi = {A/^, A/j}, where A/^ and \j.^ are given 
as in (3.8). Then we have known that Hf^p^f = 1 and that there exists a right-coprime 
factorization of the plant over Af for every / G Fi. By Lemma 4.2, the transposed plant P* 
is stabilizable. We can construct its stabilizing controller by analogy to Example 5.2. Further 
we see that for both A/j and A/j, the transposed plant P* has right-coprime factorizations 
over Axj and Axj ; that is, P has left-coprime factorizations over Axj and Axj^ . Thus let 
F2 = {A/i, A/2}. For A/j G F2, we have the matrices Nj-^ Di^, F/^, X/^ over Axj such that 
Nj.Yj, + Nj.Xj, = E2 and 



Ni, = [ll Y:, = [l 0], X,, = [l\] 



D 



h 



Xj^^ai,{{l + z){l + 2z){l-3z)) 

Xj_^ai,{{l + z){l-3z){l + 2z + 4z^)) 1 
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On the other hand, for A/j G F2, we have the matrices Ni^ Di^, Yj^, Xj^ over Axj,^ such that 
Ni^Yi^ + Ni^Xi^ = E2 and 



Nj, = [ll ¥,, = [! 0], X,, = [l\] 



D 



h 



\J^^aj,{{l + z){l + 2z){l-3z + z^)) 

1 -A7>/,((1 + 2z){l + Z + z^){l -3z + z^)) _ 



Now we let F = {Xl^hh,, AfJ (= {fj^ \ f^ ^ f,,^ e F^}). Then F still generates A 
since A?^ +2\i^\i^ + \% = 1. 

In the following we consider the case f = Xj . Then using the matrix K of (4.3), we 
have (5.1) with u = 1, Kf = K and V} = A/.A/.T. 

Then the ideal ImAif^f) i^ generated by 

Xl, aj,\l{il + z){l-3z){l + 2z + Az'')), aj.XlHl + z){l + 2z){l - 3z)). 

Thus since each of them is in Ap/^, (5.2) holds. Further we can observe that for any integer ^ 
greater than 1, (5.3) holds since Af G ImAifKf). 

For the other cases / = A/^ A/2 and f = Xj^, we can follow the relations of (5.2) and (5.3) 
analogously. Details are left to interested readers. D 

Remark 5.4. Since Anantharam's example in [1] is artificial, we do not present here the 
construction of a stabilizing controller. However we can construct it as part (c) in the proof 
of Theorem 3.3 (Since Anantharam in [1] did not consider the causality, we let Z = {0} so 
that V = J^). 

6. Causality of Stabilizing Controllers. In this section, we present two facts: (i) for 
a stabilizable causal plant, there exists at least one stabilizing causal controller and (ii) the 
stabilizing controller of the strictly causal plant is causal, which inherits Theorem 4.1 in § III 
of [14, p.888] and Proposition 1 of [11]. 
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Proposition 6.1. For every stabilizable causal plant, there exists at least one stabiliz- 
ing causal controller of the plant. 

Proof. In the construction of the stabilizing controller in part (c) of the proof of Theo- 
rem 3.3, the denominator matrix of (5.6) is Z-nonsingular. Suppose that the obtained sta- 
bilizing controller is expressed as B^^A with the matrices A and B over A such that B 
is Z-nonsingular. Then since the relation B^^A = {det{B)Em)^^ia.dj(B)A) holds, every 
entry of B'^A is causal. D 

Proposition 6.2. For every stabilizable strictly causal plant, all stabilizing controllers 
of the plant must be causal. 

Proof. Suppose that the plant P is stabilizable and strictly causal. Suppose further that C 
is a stabilizing controller of P. We employ the notation from part (c) of the proof of Theo- 
rem 3.3. Thus, a/„, A/„ G A\Z and F/^iV/^ + X^„Z)/„ = E^ with P = Ni^^Dj^^ e -F(^a,J 
from (5.5). Let Zx^ = {z/l ■ u E Axj \ z E Z,u is a unit of Axj }. Then this Zxj is again 
a principal ideal of Axj . 

Observe here that Lemma 8.3.2 of [13] and its proof hold even over a general commuta- 
tive ring. According to its proof, there exist matrices A and B over Axj such that C = B'^A 
and ANJ^^ + BDj^ = E^ (A and B correspond to T and S, respectively, in the proof of 
Lemma 8.3.2 of [13]). Observe also that every entry of Ni^^ is in Zxj . Thus reviewing the 
proof of Lemma 3.5 of [14], in which the calligraphic H and K in [14] correspond to Axj 
and Zxj , respectively, we have det B E Axj \Zxi ■ This implies that B^^A E p™^" by 
noting that A/(, E A\Z. Thus C is causal. Q 

7. Further Work. In this paper we have presented criteria for feedback stabilizability. 
We have also presented a construction of a stabilizing controller to which Sule's method can- 
not be applied. Recently the first author[7] has developed a parameterization of stabilizing 
controllers, which is based on the results of this paper and which does not require coprime 
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factorizability. This can be applied to models to which Youla-Kucera parameterization can- 
not be applied. 

REFERENCES 



[1] V. Anantharam, On stabilization and the existence of coprime factorizations, IEEE Trans. Automat. 
Control, AC-30 (1985), pp. 1030-1031. 

[2] M. F. Atiyah and I. G. MacDonald, Introduction to Commutative Algebra, Addison-Wesley, Read- 
ing, MA, 1969. 

[3] C. A. Desoer, R. W. Liu, J. Murray, and R. Saeks, Feedback system design: The fractional 
representation approach to analysis and synthesis, IEEE Trans. Automat. Control, AC-25 (1980), 
pp. 399-412. 

[4] K. O. Geddes, S. R. Czapor, and G. Labahn, Algorithms for Computer Algebra, Kluwer Aca- 
demic Publishers, 1992. 

[5] J. D. Greenfield, ed.. Microprocessor Handbook, John Wiley & Sons, 1985. 

[6] B. R. McDonald, Linear Algebra over Commutative Rings, vol. 87 of Monographs and Textbooks in 
Pure and Applied Mathematics, Marcel Dekker, 1984. 

[7] K. Mori, Parameterization of stabilizing controllers over commutative rings, 1999. (Submitted to SIAM 
J. Control and Optimization). 

[8] K. Mori and K. Abe, Feedback stabilization and robustness of stabilizability over integral domains, 
inMTNS'96, 1996. 

[9] , Improvement of generalized elementary factors and new criteria of the feedback stabilizability 

over integral domain, in European Control Conference ECC '97, TH-M D3 (533), 1997. 

[10] S. Shankar and V. R. Sule, Algebraic geometric aspects of feedback stabilization, SIAM J. Control 
and Optimization, 30 (1992), pp. 1 1-30. 

[11] V. R. SULE, Feedback stabilization over commutative rings: The matrix case, SIAM J. Control and 
Optimization, 32 (1994), pp. 1675-1695. 



39 



[12] , Corrigendum: Feedback stabilization over commutative rings: The matrix case, SIAM J. Control 

and Optimization, 36 (1998), pp. 2194-2195. 

[13] M. ViDYASAGAR, Control System Synthesis: A Factorization Approach, MIT Press, Cambridge, MA, 
1985. 

[14] M. ViDYASAGAR, H. SCHNEIDER, AND B. A. FRANCIS, Algebraic and topological aspects of feedback 
stabilization, IEEE Trans. Automat. Control, AC-27 (1982), pp. 880-894. 



